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Abstract. Recent measurements of type la supernovae 
indicate that distant supernovae are substantially fainter 
than expected from the standard flat cold dark matter 
model. One possible explanation is that the energy density 
in our universe is in fact dominated by a cosmological 
constant. Another possible solution is that there are large 
amounts of grey dust in the intergalactic medium. Dust 
grains can be grey either because they arc non-spherical or 
very large. We have numerically investigated whether grey 
dust can be emitted from high redshift galaxies without 
also emitting standard, reddening dust, which would have 
been visible in the spectra of high redshift objects. Our 
finding is that grain velocities are almost independent of 
ellipticity so that if greyness are due to the grains being 
elongated, it will not be possible to separate grey dust 
from ordinary dust. We also find that velocities are fairly 
independent of grain size, but we cannot rule out possible 
sputtering of small grains, so that large, grey dust grains 
could be preferentially emitted. Therefore, our conclusion 
is that grey dust is an unlikely explanation of the data, 
but we cannot rule it out if the grey dust consists of large, 
spherical grains. 

Key words: ISM: dust, extinction - Galaxies: intergalac- 
tic medium - Cosmology: dark matter 



1. Introduction 

In recent years the prospect of using type la supernovae 
for distance measurements in the universe has attracted a 
great deal of attention. Many studies indicate that these 
supernovae are very close to being standard candles, so 
that measuring their effective magnitude amounts to di- 
rectly estimating their luminosity distance. Two large 
scale surveys of high redshift supernovae are currently in 
progress (Perlmutter et al. 1998, 1997; Garnavich et al. 
199S| ; Riess et al. |1998| ; Schmidt et al. |1998fl , and both 



groups have published results for the determination of 
the deceleration parameter, q . Surprisingly, the results 
are statistically incompatible with the value q — 1/2 ex- 
pected in the standard, flat matter dominated Friedmann 



model (Kolb & Turner 199C ). The observed high red- 
shift supernovae are significantly fainter than expected for 
qo = 1/2. Rather they are compatible with a medium den- 
sity, flat universe, where roughly 70% of the en ergy d ensity 
comes from vacuum energy (Perlmutter et al. 1998 ) 

If this turns out to be true it is an astonishing result 
because it means that, taken together with the recent re- 
sults on neutrino masses from Super-Kamiokande (Fukuda 
et al. 1998), there are at least four components contribut- 
ing almost equally to the cosmic energy density budget: 
baryons, neutrinos, cold dark matter and vacuum energy. 
Even without vacuum energy this is a remarkable fine tun- 
ing, but including vacuum energy as a fourth component 
makes the conceptual problem substantially worse. 

For this reason it is important to check all alternative 
explanations of the supernova data that can be thought 
of. One such possibility is evolution effects in the super- 
nova progenitors. However, the data seem to indicate that 
the supernova population at high redshift is remarkably 
similar to the one known at low redshi ft (Per lmutter et al. 



1998; Hoflich, Wheeler & Thielemann 1998). Another in- 



triguing possibility which has been advocated by Aguirre 
(Aguirre 1999a , 1999b| ) is that there could be unseen dust 
contamination in the data. In the data reduction it has 
been assumed that any possible dust follows the standard 
extinction law known to govern dust in the Milky Way 
(Cardelli, Clayton & Mathis |1989| ). 

Based on this, almost no dust is found in the line of 
sight to any of the supernovae. However, this conclusion 
is based on the fact that dust in the Milky Way red- 
dens significantly, simply because it is primarily made 



up of small, spherical grains (Draine & Lee 1984; Aguirre 



1999b). This need not be the case. If there is some popu- 



lation of very large and/or elongated dust grains which is 
homogeneously distributed in the universe, this would re- 
main undetected in reddening surveys because of its very 



flat extinction curve (Aguirre 1999a, 1999b), i.e. it would 
be "grey" dust. In fact the idea of grey dust is quite 
old and has been advocated as a means of thermalising 
the cosmic microwave background radiatio n in t he steady 
state cosmology (see for instance Peebles 1993 for a re- 
view) . 



2 



Jakob T. Simonsen & Steen Hannestad: Can Dust Segregation Mimic a Cosmological Constant? 



The above scenario could perhaps be a viable expla- 
nation of the data, given that enough such dust can be 
expelled from galaxies at high rcdshift without simultane- 
ously expelling small dust grains. 

The purpose of the present paper is to quantify this 
by numerical modelling of dust expulsion. We have wished 
to investigate whether it is in fact possible to expel only 
very elliptical and large grains. We have not worried about 
the possible production mechanisms for dust grains, only 
about their dynamical behaviour. A completely different 
problem is whether it is at all possible to produce the re- 
quired amounts of large dust grains in the first place. How- 
ever, by using our relatively simple model we are able to 
make the very robust prediction that dust segregation at 
the level needed to explain the data is not possible, unless 
very special conditions prevail in the host galaxies. Note 
that in the present paper we only consider radiation pres- 
sure as a source of dust expulsion. There could be several 
other ways of removing dust from galaxies, such as winds 
from supernovae (Such kov et al. 1994 ) or galaxy collisions 
(see e.g. Gnedin 1998| ). However, these methods seem un- 
likely to preferentially emit large dust grains, and for that 
reason we have, in accordance with Aguirre ( 1999a , 1999b ) 
concentrated on radiation pressure, even though galaxy 
collisions might be the most efficient method of dust re- 
moval at high redshift (Gnedin 1998). 

The paper has been sectioned as follows: Section 2 con- 
tains a discussion of the equation of motion for dust grains 
in a galactic environment, section 3 contains our main nu- 
merical results, and finally section 4 is devoted to a dis- 
cussion of our findings. 



2. Equation of motion for dust grains in a galactic 
medium 

At present the supernovae with the highest measured red- 
shifts are at z 



0.8 (Perlmutter et al. |1998[ ). This means 
that if dust is to explain the data it must be expelled and 
distributed uniformly prior to this epoch. For this to hap- 
pen it must have been expelled from the host galaxies at 
very high redshift, probably z ~ 3 or more (Aguirre 1999a, 



1999b). Relatively little is known about the structure of 



galaxies at such distances, except that they are richer in 
gas content than our own galaxy at present. 

We have taken a heuristic approach to our lack of 
knowledge about the nature of the host galaxies and cal- 
culated dust expulsion for a range of different gas contents 
and emission temperatures. We find that our results de- 
pend very little on these assumptions. 

We assume that a given spiral galaxy can be reasonably 
well approximated by three components: disk, bulge and 
halo. The disk is taken to be an infinitely thin axially 
symmetric exponential distribution of matter. The bulge, 
which is massive and luminous, is just added to the center 
of the galactic disk, and the halo is a spherical symmetric 
mass structure considered to be completely dark. 



The three main contributions to the force acting on a 
dust grain moving through a galaxy are gravity, radiative 
forces and drag forces. As the galaxy is assumed to be 
axially symmetric, all calculations are done in cylindrical 
coordinates (i?, z, <j>). 

In the following we shall always work with ellipsoidal 
dust grains, with semiaxes a and b. The effective radius of 
such a grain is then defined by the relation 



4 o 4 2 



(1) 



2.1. Gravitational forces 



The gravitational force on a dust grain of mass m g = p g V, 
where p g is the grain density, is simply given by 



F G = m g G(r), 



(2) 



where G(r) is the gravitational field intensity at the point 
r. 

The disk has an exponential mass distribution with 
surface density 



E(iZ) = E exp(-i?/i? d ), 



(3) 



where Rd is the disk scale length. The potential generated 
by this disk is most easily expressed via Bessel functions, 



and takes the form (see e.g. Binney & Tremaine 1987) 



<b(u r\- vrv u 2 r MkR)cM-kz) 
*(R,*)--2*G2°RiJ Q [1 + (fci?d)2]3 /2 



dfc, 



(4) 



J n (x) being the Bessel function of the first kind of order 
n. The gravitational force from the disk on a dust grain 
with mass m g is thus given by 

p2 [°° Ui(kR)exp(-kz) j U 
-Fgrav.disk.fl = -m s 2nGT, R d / — _ ,„,„ ,„ dfc, 



F, 



grav,disk,z 



o [1 + (kR d ) 2 } 3 / 2 
2 [°° fcJ (fcE)exp(-fcz) 

~ ms27rGS0 M [1 + (*W 2 



dfc. 
(5) 



The potential from the bulge is that of a point mass 
with mass Mb. That is, 

GMbTOg 



F 



grav, bulge, r 



(6) 



where r = ^ R 2 + z 2 . Projecting on to the R and z axis 
gives 



F 



GM h m„ R 



R 



grav, bulge,/?. 



- ■ ~ - GM bm g(i?2 + ^ 2)3/2 . 



Fc 



rav, bulge, z 



r 

GM^rricr z 



— = —GMbm g 



(R 2 + z 2 ) 3 / 2 ' 



(7) 



For the halo we assume an isothermal sphere distribu- 



tion 



ai + r z 



(8) 
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Integrating the equation for p yields 

M(r) = inpoa^lr — ah arctan(r/a n )], (9) 

where po is the central density and ah is the halo mass 
scale length. Setting Mh = M(r), this results in a central 
force, which after projecting on to the R and z axis takes 
the form 

R 

-Fgrav,halo,ii = — GM^nig 



-Fgrav,halo,z = — GM^m^ 



' (R 2 + Z 2 ) 3 / 2 ' 

z 

(R 2 + z 2 ) 3 / 2 ' 



(10) 



2.2. Radiative forces 



The radiation pressure force can be written in a way com- 
pletely analogous to the gravitational force (Ferrara et al. 



1991) 



First, we shall assume that the spectral function, VL V , 
of the emitted light is constant throughout the galaxy, i.e. 



n u (R,z) = si v . 



(ii) 



If this is the case then we can calculate a spectrally aver- 
aged radiation pressure coefficient for a given dust grain 
as 



dvQ. v Q pT (v), 



(12) 



which applies at all positions in the galaxy. In practise we 
shall assume that the galaxy emits light as a black body 
with some temperature Tq, so that Q. v = _B„(Xo), where 
B v is the Planck function, and Q* r (a c fr) = Q P r(aeff ; 2o). 

The disk luminosity distribution is assumed exponen- 
tial like the mass distribution, and with the same scale 
length i?d 



I(R) =I exp(-R/R d ). 



(13) 



From these simple assumptions the radiation pressure 
force from the disk is given as 

2ir 2 a 2 cS Q; r J f°° kJ x (kR) exp(-fe) 



ad, disk, R 



F 



2ir 2 a 2 cS Q*h 



rad.disk.z 



[l + (fci? d ) 2 ] 3 / 2 

kJ (kR) exp(— /cz 
[1 + (fci? d ) 2 ] 3 / 2 



-dk, 



■dk, 



analogous of Eq. (||) . 

The force from the bulge is similarly given as 



ad,buls;c,r 



7Fa cff ( 3pr ^ 



(14) 



(15) 



c Anr 2 

where Lb is the total luminosity of the bulge. In the R and 
z directions this is 

_ "cffQpr^b 
-^rad, bulge, R — 



R 



4c 



(r 2 + z 2 yi 2 ' 



Jrad, bulge, z 



C[ eff < 3pr^'b 



4c 



(R 2 



2)3/2- 



(16) 



again, apart from numerical factors equal to Eq. (Q). 

Note that there is an important issue which we do 
not address here, namely the opacity in the disk itself 
(Davies et al. 1998). Since the disk contains significant 



amounts of dust it will obscure the starlight and lead to 
smaller radiation pressure. This was determined to be an 
important effect for the dynamical evolution of grains by 
Davies et al. 



(199S) 



However, it is not entirely clear how disk opacity affects 
the grain evolution. It leads to smaller radiation pressure, 
but on the other hand it also reddens the light, leading 
to a lower effective temperature of the radiation field. We 
model this behaviour by calculating grain evolution for 
several different emission temperatures. 

2.3. Rotation 

The stars and gas in the disk of spiral galaxies have purely 
rotational orbits. This means that the rotational velocity 
is given by 



Vrot 



R 



OR' 



(17) 



Initially the grain also rotates with this velocity, but as the 
grain moves outwards the rotational velocity decreases. 
If gas drag is neglected the angular momentum in the z 
direction is conserved so that 



v m t(R) = v m t(Ri)(Ri/R). 



(18) 



However, for typical spiral galaxies the rotational velocity 
of the stars remains constant out to very large values of 
R. This means that grains will, in general, be spun up 
by interaction with the disk gas. This effect is in practice 
quite small, but we have included it anyway. 

2.4- Drag forces 

The interaction of dust grains with gas plays an important 
role in grain evolution. Assuming that the grain moves fast 
compared with the thermal velocities in the gas, the exact 
expression for the drag force experienced by a grain is 
(Il'in |1994|) 



drag 



= -Tiabp g&s v 2 —E([l - e 2 ] sin 2 9)v g 



(19) 



where a and b are the semiaxes of the spheroidal grain, 
e = a/b, v g the grain velocity relative to the gas, E(x) 
the complete elliptic integral of the second kind and the 
angle between magnetic field direction and wave vector of 
radiation. 

Interstellar dust grains rotate with large angular ve- 
locity because of the angular momentum transferred to 
them by collisions. Non-spherical grains are likely to ro- 
tate along the axis of maximum moment of inertia, with 
the angular momentum being aligned along the magnetic 
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field lines (Il'in 1994). For this reason O is likely to be 
small so that the argument to the elliptic integral is close 
to 0. Since -E'(O) = 7r/2 the drag force becomes 



-nabp 



gasfgVg. 



(20) 



As mentioned above, this drag force also causes the grain 
to be spun up by the gas, so that, as long as it stays in 
the gas, , is driven towards the rotational velocity of the 
gas. Note that [1-e 2 ] sin 2 6 G [0, 1] so that E{x) <E [1, tt/2] 
always, i.e. choosing x = produces at worst an error of 
50% in the drag force. 

Also, the above formula only applies if the grain veloc- 
ity is larger than the thermal velocity of the gas. However, 
using Eq.([l9|) leads to small errors, as will be discussed in 
section 3.1. 



2.5. Sputtering 

In general, grains are eroded as they pass through the 
galactic medium, by collisions with gas particles a nd by 
radiation (Drai ne fc Salpeter 1979 , Ferrara et al. 1991 , 
Shustov & Vibe 1995 ). The amount of sputtering increases 



with gas density, both because of the higher collision rate 
and because the grain spends more time in the region with 
high gas density. The effect of sputtering is primarily to 
destroy very small dust grains, but leave large dust grains 
more or less unscathed (Ferrara et al. 199ll Davies et al. 



1998| ). It is very difficult to assess the amount of sputtering 



in a given galaxy. For galaxies similar to the Milky Way, 
Ferrara et al. (1991) find that sputtering is not important 
unless the grains are very small (a ~ 10 — 20 nm), whereas 
Shustov & Vibe (1995) find that grains up to 50 nm could 
be destroyed. We can therefore not rule out possible de- 
struction of grains up to about this size. 

2. 6. Magnetic forces 

As has been noted many times in the literature, galactic 
magnetic fields are important for the dynamical evolution 
of dust. However, very little is known about the nature of 
these fields, so we have chosen to neglect possible magnetic 
fields in the present study. In general, the effect of the 
magnetic field will be to enhance the confinement of grains 
because the drag force on grains is enhanced (Ferrara et 
al. 1991) and because they are forced to move along the 
field lines (Ferrara et al. 1991; Davies et al. 1998[ ). 



3. Numerical results 

In order to do actual numerical calculations it is neces- 
sary to specify relevant parameters for the host galaxy. In 
practice we have taken the present day Milky Way with 
parameters as specified in Table 1 (Binney & Tremaine 



1987, Ferrara et al. 1991). This may not be a very accu- 
rate description of a typical high redshift galaxy, but as 
will be shown the numerical results are fairly insensitive 



Table 1. Physical parameters for our sample galaxy (the 
Milky Way). 





Parameter 


value 


Disk 


So 


381.9 M pc"^ 






5 kpc 




h 


108.28 L pc" 2 


Bulge 


M b 


1.0 x 10 10 M Q 




L b 


1.5 x 10 9 L Q 


Halo 


Pa 


6.3 x 10~ 24 g cm -3 
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Fig. 1. Grain velocity v and displacement r for 
the present day Milky Way model, with (a e g, e) = 
(50 nm, 1), (50 nm, 1/8), (300 nm, 1) (full, dashed and dot- 
ted curves respectively) , and initial galactocentric distance 
Rq = 2 kpc (upper three curves) and Rq = 12 kpc (lower 
three curves). Notice that the initial velocity is due to 
rotation. 



to these assumptions. The most important issues are the 
emission temperature of the galaxy and the amount of gas 
contributing to the drag force. We have modelled the radi- 
ation field of the galaxy as a black body with temperature, 
T = 6000 K, 9000 K or 12000 K. 
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Fig. 2. Grain velocity v and displacement r for 
our extremely gas rich model, with (a e ff j e) = 
(50 nm, 1), (50 nm, 1/8), (300 nm, 1) (full, dashed and dot- 
ted curves respectively), and initial galactocentric distance 
i?o = 2 kpc (lowest initial v) and Ro = 12 kpc (highest ini- 
tial v). Notice that two of the Ro = 12 kpc curves overlap, 
and that the initial velocity is due to rotation. 

Furthermore we have considered three different models 
for the gas distribution. In all three cases we model the gas 
by two different components, a cold gas and a hot halo gas. 
Both components are assumed to have the number density 
distribution 



n(R, z) 



n e 



zjz 



(21) 



i.e. constant density as a function of R and an exponen- 
tial fall-off with z. The cold component is assumed to have 
z s — 170 pc, and the hot z s — 2 kpc. For the three differ- 
ent models we then use the following values for no 



^o,hot = 1-0 x 10~ 4 cm~ 3 , n 0iCO id = 2.0 x 10~ 3 cm -3 
^o,hot = 1-0 x 10~ 2 cm~ 3 , n 0i coid = 2.0 x 10 _1 cm -3 
n ,hot = l-° x 10_1 cm -3 , n ,coid = 2.0 cm~ 3 
The first model corresponds to the extreme case of an al- 
most empty galaxy, the second roughly to the present day 



Milky Way (Lockman \L99l\ Heiles |l99l| ), and the third to 
an extremely gas enriched galaxy. In all cases we assume 
that gas is made entirely of hydrogen. Note that we have 
assumed a model with constant density as a function of R. 
This may seem a strong assumption, but as will be shown 
in the next section the numerical results depend relatively 
little on it. 

The next important iss ue is th e physical composition 
of the dust grains. Aguirre 1999a assumed pure graphite 
grains, b ut rela xed that assumption in a subsequent study 
(Aguirre 1999b ), where a more realistic galactic dust com- 
position of silicate and graphite was assumed. 

Again, we have taken an heuristic approach and calcu- 
lated our models for two different cases, amorphous car- 
bon and graphite-like carbon; silicates are expected to be- 
have very similarly. To find the radiation pressure coef- 
ficient for t hese g rains we have used dielectric data from 
Jager et al. 1998. They have obtained dielectric data on 



carbon pyrolised in the laboratory at different temper- 
atures. Their 400 C sample has an amorphous structure 
while their 1000 C sample closely resembles graphite. Both 
are bulk samples and therefore isotropic. Note that the 



graphite needles used in the Aguirre ( 1999a ) model are 
highly anisotropic, i.e. they are macrocrystals, whereas 
our data comes from a bulk sample containing microcrys- 
tals. However, in practice the difference between these two 
assumptions is likely to be small, so that our results ap- 
ply to the Aguirre (1999a) scenario. We take these two 
samples as extreme cases of carbon composition in dust 
grains. Using this data we have proceeded to calculate the 
radiation pressure coefficient by use of the Discrete Dipole 
approximation as implemented in the DDSCAT package 
(Draine & Flatau 1998). The resulting Q* T are tabulated 
in|£J 

Next, we have performed the full numerical solution 
to the dynamical equations for ellipsoidal grains of dif- 
ferent ellipticities and sizes. Specifically we have calcu- 
lated for a cB = 25, 50, 100, 150, 200, 250, 300 nm, and e = 
1,1/2,1/4,1/8,1/16 for each a c g. The calculations have 
been done for twelve different initial values of R, equally 
spaced from 1-12 kpc, and always starting at z = 0.2 kpc. 
The grains were followed for a period of 2 Gyr to deter- 
mine their final velocity. Figs, [l] and || show some examples 
of velocities and displacements r — (R 2 + z 2 ) 1 / 2 as a func- 
tion of time, for the present day Milky Way and for our ex- 
tremely gas enriched model. In each case the velocity and 
displacement are shown for three combinations of a e g and 
e, namely (a e g,e) = (50nm, 1), (50nm, 1/8), (300nm, 1), 
and for two different initial galactocentric distances Ro — 
2, 12 kpc. We have used Q values for the appropriate size 
and shape at To = 9000 K. Both figures show that it is 
quite difficult to expel dust which has a large initial galac- 
tocentric distance, whereas most dust is expelled rapidly 
from the galactic center. In the gas rich model, we find 
that grains starting at R — 12 kpc are almost completely 
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Fig. 3. Con tours of constant velocity in the (a e g, e ) plane for graphite like carbon grains (1000 C sample from Jager 
et al. |1998[). 



confined by drag forces. The figures also show that grains 
of different sizes and ellipticities behave fairly similarly. 

After solving the dynamical equations, we have then 
calculated the final velocity of the dust grains, averaged 
over the entire galaxy. To do this we have assumed that 
the dust production rate is proportional to the light profile 
of the disk, i.e. R cx e~ R / Rd . In some cases dust grains are 
confined within the galaxy by drag forces in the gas. This 
happens particularly in the outer parts of the host galaxy 
and in the gas rich model. For such grains we set the 
final velocity equal to zero. This means that our averaged 
velocity is in fact an average over all grains, not only the 
ones which escape. Thus, the averaged final velocities are 
in fact a folding of the velocity of escaping dust grains and 
escape probabilities. 

Figs. U and [I] summarise our results for graphite and 
amorphous carbon respectively, showing contours of con- 
stant velocity in the (a e s, e _1 ) plane, for all combinations 
of base density no (represented by no, hot on top of the 
columns) and temperature To (to the right of each row). 
Exactly as expected the averaged velocity is a decreasing 



function of the gas density and an increasing function of 
emission temperature (unless the grains are very large). 

Our main result is that the final grain velocity is a very 
slowly varying function of ellipticity. The reason is that al- 
though the effective Q-value increases with ellipticity if the 
effective radius is kept fixed, the drag force increases too. 
Drag therefore acts as a "velocity-moderator" in the sense 
that it tends to erase velocity difference due to ellipticity. 

The differences in final velocity due to varying grain 
size is straightforward to understand. The acceleration due 
to radiation pressure is proportional to a~g so that one 
should expect decreasing final velocities for larger grains. 
However, the effective Q- value is a strongly increasing 
function of a c g for small grains (up to roughly 100 nm), so 
that in the end there will be some a e s where the velocity 
peaks. Exactly this effect is seen in Figs. 3 and 4. 

Note that in the upper left plot of fig. ^ the veloc- 
ity increases again at very small radii. This is probably a 
numerical artifact produced by uncertainty in the calcu- 
lation of the effective Q-value by the program DDSCAT. 
A few test runs have shown that increasing the number 
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Fig. 4. Con tours of constant velocity in the (a G ff,e 1 ) plane for amorphous carbon dust grains (400 C sample from 
Jager et al. |l998j). 



of discrete dipoles in the grain does diminish the Q value, 
but at much incresed computational time. 

Our numerical results thus indicate that it will be dif- 
ficult to separate dust grains of different ellipticities, but 
that in certain cercumstances it could be possible to expel 
large grains only. 

3.1. Model dependencies 

We have investigated how dependent our conclusions are 
on our assumptions about the drag force and the number 
density distribution of the gas, by performing some sample 
calculations using more conservative assumptions. 

In calculating the drag force, the grain velocities were 
assumed to be large compared to the thermal velocities 
in the gas. But adopting gas temperatures of 100K and 
10 6 K for the cold and hot gas component respectively, 
the average thermal velocity of the hot component is 
u gas ~ 130km/s, which is not negligible compared to the 
smaller grain velocities. In order to check the validity of 
using our Eq. ( |l9|) we recalculated the final averaged ve- 



locities of amorphous carbon grains in the model with base 
densities representative of the present day Milky Way (the 
second model) , using the approximate formula for the drag 
force given by Il'in (1994, Eq. 23). This approximation is 
designed to be valid for all grain velocities. The results are 
shown in the first column of figure ||. 

Even though our equation can produce errors of 
roughly 25% in absolute velocities, it has no effect on the 
conclusion that grain velocities are fairly independent of 
size and ellipticity. 

For the gas we assumed that the number density dis- 
tribution was constant as a function of R. We have inves- 
tigated the effect of this by performing a calculation using 
an exponential fall-off with both R and z. We use base 
densities of no iCO id = 1.1 cm" 3 and rio,hot = 0.05 cm" 3 , 
and the same scale length as the stellar disk. This pro- 
duces densities equal to those of the present day Milky 
Way at Rq ~ 8.5 kpc. Here we return to using the expres- 



sion for the drag force given by Eq. (] 
shown in the second column of figure [ 



3). The results are 
Again, this yields 



results which are somewhat different in absolute magni- 
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Fig. 5. Contours of constant velocity in the (a c ff,e 1 ) plane for amorphous carbon dust grains using the approximate 
drag force from Il'in (1994) (left column), and with exponential fall-off in number density with R (right column). 



tude relative to those found from the constant density gas 
disk. However, the velocity differences due to different size 
and ellipticity are still small, which is the main point. 

In light of the large uncertanties in these calculations, 
because of our poor knowledge of the nature of high- 
redshift galaxies, we think our approximations are rea- 
sonable. They do yield different absolute results, but have 
no bearing on our conclusion about the relative velocities 
of different dust grains. 



4. Discussion 

We have performed numerical calculations of how dust 
grains are expelled from galaxies by the radiation pres- 
sure produced by starlight. The calculations were done 
for a whole range of different sizes and shapes, as well as 
for several different host galaxy types. In all cases, final 
grain velocities are almost independent of ellipticity. Our 
conclus ion therefore is that the original scenario of Aguirre 
(|l999a|) does not work. 



Also, if the galaxy does not contain significantly more 
gas than the present day Milky Way, the relative velocity 
difference between grains of different sizes is quite small, 
unless the grains become small. 



The mechanism proposed by Aguirre (1999t) requires 
that only dust with effective radius larger than about 
100 nm should escape to the intergalactic medium, in or- 
der not to have excessive reddening (this can be relaxed 
to about 50 nm if only graphite grains are considered 
(Aguirre 1999b)). Our results indicate that this mecha- 
nism does not work for the two models with the least gas 
content (the first two columns in Figs. || and ^), unless 
the smaller grains are destroyed completely by sputtering. 
However, the findings of Ferrara et al. ( |199l[ ) and Davies 
et al. (1998), indicate that sputtering is not effective un- 
less the radius becomes as low as 10 nm. On the other 



hand, Shustov & Vibe (1995) find that grains up to 50 



nm can be destroyed. In light of our poor understanding 
of the sputtering mechanism it is not possible to rule out 
the Aguirre model. However, just from looking at grain 
velocities the model does seem an unlikely explanation. 
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As for the third model with extremely large gas con- 
tent, there is a difference between small and large grains, 
the smaller grains having lower velocities. This seems to 
be exactly the type of selection effect we are looking for 
and, furthermore, the model where it happens could be 
quite similar to high redshift galaxies. However, the result 
should be treated with caution because this model is al- 
ready close to the limit where all grains are completely 
confined by gas drag. Altogether we find that there are 
models where dust segregation can be produced, but that 
the physical conditions have t o be jus t right. We can there- 
fore not rule out the Aguirre ( 1999b ) model, only say that 
it requires some fine tuning. 

Our results indicate that most of the dust produced in 
galaxies will be expelled by radiation pressure. Of course 
we have neglected galactic sputtering of grains, which 
might destroy very small grains. However, measurements 
indicate that the intergalactic mediu m in fact has a fairly 
high metal content (Mushotsky 1996). This in turn is likely 
to mean that most of the produced dust is expelled from 
the host galaxies. However, it also seems very likely that 
almost all of the dust is subsequently destroyed in the in- 
tracluster medium (Dwek et al. 1990, Stickel et al. 1998| ). 
Note that since there is probably effective sputtering in 
the intracluster medium, this could also mean that selec- 
tion of large dust grains could be due to sputtering in the 
ICM. This also limits our ability to rule out the Aguirre 
( |l999b[ ) scenario. We can only say that sputtering within 
galaxies is an unlikely explanation. 

The above discussion does not take into account that 
radiation pressure may not even be a dominant mechanism 
for dust expul sion at high redshift. As mentioned above, 
Gnedin ( 1998 ) has argued that galaxy mergers may pro- 
vide the most efficient mechanism for expulsion of grains 
into the intergalactic medium. If this is the case there is 
even less reason for suspecting that grains with flat opac- 
ity curves should preferentially be emitted, and the above 
scenario will be even less credible. 

Note added - After this paper had been submitted 
there has been a growing interest in the possibility of pro- 
genitor evolution systematically biasing the SNIa data. 
At present it seems unclear what the effects of evolution 
are, but it has been pointed out that if there is evolution, 
then it will be very difficult to extract cosmologic al pa- 
rameters from the data (Drell, Loredo & Wasserman 
Dominguez et al 



Table A.l. Effective radiation pressure coefficients, Q* r , 
for different grain sizes and ellipticities. 

Carbon, T = 6000 K 



1999 



Riess et al. 1999; Wang 1999) 



1999 
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Appendix A: Radiation pressure coefficients 

The following table shows the averaged radiation pres- 
sure coefficients Q* r as a function of a e fi, e and radiation 



field temperature T , calculated for amorphous carbon and 
graphite. a e g is in units of lOOnm. 
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